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FACTORIZATIONS RELATED TO THE RECIPROCAL PASCAL 

MATRIX 

HELMUT PRODINGER 


Abstract. The reciprocal Pascal matrix has entries ^■ Explicit formula for 

its LU-decomposition, the LU-decomposition of its inverse, and some related matrices 
are obtained. For all results, g-analogues are also presented. 


1. Introduction 


Recently, there has been some interest in the reciprocal Pascal matrix M, dehned by 

-1 


Mij 


% + j 

3 


the indices start here for convenience with 0, 0, and the matrix is either inhnite or has 
N rows and colnmns, depending on the context. 

Richardson |6] has provided the decomposition S = GMG, where the diagonal ma¬ 
trix G has entries Gi^i = , and S is the super Catalan matrix 1210] with entries 


(2i)!(2j)! 

We want to give an alternative decomposition of M, provided by the LU-decomposition. 
We will give explicit expressions for L and U, dehned by LU = M, as well as for L~^ 
and U~^. 

Since there is also interest in M~^, in particular in the integrality of its coefficients, 
we also provide the LU-decomposition AB = M~^, and give expressions for A, B, A~^ 
and B~^. 

In the last section, we provide g-analogues of these results. 


2. Identities 

The LU-decomposition M = LU is given by 


and 


^ {-iyj\j\i\{i - 1)! 


for i > 1. 


For z = 0, the formula is Uqj = 1. 
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The formula that needs to be proved is 


0 <A:<min{ 2 ,ji} 


i+j 
3 


-1 


which is equivalent to 

2i\i\j\i\ 


1 + 




(2t)!(2j)! 


l<k<mm{ij} 

The von Szily identity [3 |21 E] is 

mm 

i\i\{i + j)\ 


i + j 
3 


-1 




k£Z 


2i 


2j 


i + kj\j + kj' 


and an equivalent form is, by symmetry, 


il3l{i + j)\ 


k>l 



Thus, the identity to be proven is now 

i + j'\ + j)l 

3 J m-im + 

which is obviously correct. 

The formula for L~^ is for i > j > 0: 

1 ^ (-l)^-^'z!z!(i + i - 1)! 

+ (2^ -1)^(^ - jy.jyi • 

If necessary {i = j = 0), this must be interpreted as a limit. 
To check this, we consider 

iH\{2k)\ {-imk\k\{k + j -1)\ 
z + k)\{i - k)\k\k\ {2k - l)!(fc - j)\i\j\ 


= 1 , 




2i\i\{-iy 


Ett- 


k 






The sum can be evaluated by computer algebra (or otherwise), and the result is indeed 
[z = j], as desired. 

The formula for U~^ is for j > z > 1 

(-l)*(j + z)!(2j)! 


and for z = 0: 


(i - + «)Ki - 

(2j)! 


Uo} = 


3'-3'- 


The fact that J^k^hk^kj = [* = j] can also be done by computer algebra. Since 
there are a few cases to be distinguished, it is omitted here. 
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The LU-decomposition AB = M ^ depends on the dimension N and is given by 

_ -j- l)!j!(A^ + i - 1)! 

(-l)^W(iV + j-l)! 

Since M~^ does not have “nice” entries, we rather provide formulae for A~^ and B~ 
and prove the identity B~^A~^ = M instead. The results are: 

“ (j-f)!(A^ + i-l)! ■ 

First we prove that these are indeed the inverses. We consider 

(_l)*-fc(Al - k - l)\kl{N + i - 1)\ {N - j - l)\j\{N + k - 1)\ 


E 


t\{N - i - 1)\{N + k - l)\{i - k)\ k\{N - k - 1)\{N + j - l)\{k - j)\ 
fN + z-l)\{N-j-l)\j\ ^ _(-1)^ 

V ) r AT T , • 1\|-| /• 


{N-i- 1)!(A^ + j- l)!d (^ - k)\{k - j)\ 


{N - i - 1)\{N + j - - i)\ 


E(-i) 


j<k<i 


i — k 




= 


{N - i - 1)\{N + j - - j)\ 

which proves AA~^ = I. Similarly 

,(-l)^+^-i(iV - 1 - i)\k\i\ {-iy+^-\N + j- 1)! 
(A:-f)!(A^ + i- 1)! j\{j - k)\{N - j - ly.kl 

AN+ j-1)\ ^ {-If 


E- 


= (-i) 


(A^ + f-l)!j!(A^-i-l)! ^ {k-i)\{j-k)\ 


E 






’ \j-k 


{N - 1 - i)W.{N + j - 1)! 


= j 1 = P = il, 


(7V + z-l)!j!(iV-j-l)!(j-*)! 

which proves B~^B = I. 

Now we compute an entry in B~^A~^: 

1 {N-j- l)!i!(iV + k-l)\ 


E- 


(A: - f)!(A^ + f - 1)! A;!(A^ - A: - 1)!(A^ + i - 1)!(A: - i)! 

(iV - 1 - ^)H!(iV - J - l)!j! ^ (-l)'=(iV + A; - 1)! 


= (-l) 


E 


(iV + i - l)!(A^ + j - 1)! ^ (A:-f)!(A^-A:-l)!(A:-i)! 



4 


HELMUT PRODINGER 


= (-l) 




{N + i- 1)! 


<■ \ / JV -I-ft. ± 

’\N-l-kJ\N-l+j 


i\j\{N - j - 1)\ ^ k i - 1 - k \ kN + k - 1\ 
{N + i- 1)! ^ Vat -1 - fcy Va^ -1 + j7 

i\j\{N — j — 1)! ki — 1 — k\kN + k — 1 


{N + i- 1)! 


E 


i-N )\N-l + j 


i\jl{N - j - 1)\ ki - 1 + N\ 

““(F+T^ryrv t+j ) 

i\j\ 

= r , -M = 

{i + j)\ 

as claimed. 

Now we use the form M~^ = AB and write the (hi) entry: 

(AT - h - l)!h!(A^ + i - 1)! {-iy+N-i(^N + i - 1)! 

i\{N - i - 1)!(A^ + k- l)!(i - k)\ j\{j - k)\{N -j - l)!fc! 


E; 


(A^ + ^ - l)!(A^ + i - 1)! 


— i — 


(A^-h-l)! 


z!(A^ - i - l)!i!(Ar - j - 1)! ^ (A^ + A; - l)!(z - k)\ {j - k)\ 


N-i'^^+j-i)^ (- 17 +^ 


J 


_^fN + i-l\l'N-k-l 
i-k )\ j-k 


0 <k<m.m{ij} 

From this representation, it is clear that this is an integer. This was a question which 
was addressed in the affirmative in [B]. 

3. ^-ANALOGUES 

In this section we present g-analogues. Dehne {q)n ■= (1 — g)(l — ... (1 — g"), and 


(?)r 


{q) k{Q)n—k 


these dehnitions are standard, see [T]. Then we have the following results. 

{Q)i{q)i{q)2j 


{q)i+j{q)i-j{q)j{q)3 




{q)i+j{q)j-i{q)2i 
-1 _ q"^"~^^^^{-^y~^{q)i{q)i{q)i+j-i 






{q) 2 i-i{q)i-j 

.l)7-i"-T+*h+i)/2(g)^.^.i(g)27g)hg)i 


{q)j-i{q)M)i-^ 


for i < i, L._. = 1, 
for i > i, 
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u. 


-1 


(_l)^g»(3» + l)/2(g)^.(g)^. 


Uo,o 


1, 


^-1 _ g(^-i)(^-^+i) {q)N-j-i{q)N+i-i{q)j 
{q)N-i-i{q)N+j-i{q)i{q)i-j 

(_1).+N-1^,0+1)/2-(N-,-1).-A(A-1)/2(^)^_.(^)^^._^ 

{q)N-^-l{qUq)^ 

Note that for g —>■ 1, we get the previous formulae. We do not discuss proofs here, 
since Zeilberger’s algorithm (aka WZ-theory) |5] proves all these results (which were 
obtained by guessing), using a computer algebra system (such as, e. g.. Maple). 


Remark. Richardson’s decomposition S = GMG still holds when all binomial 
coefficients are replaced by the corresponding Gaussian g-binomial coefficients. 
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